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Abstract 

The elastic scattering of the 16 0+ 28 Si system has been analyzed with a mod- 
ified potential within the framework of the optical model over a wide energy 
range in the laboratory system from 29.0 to 142.5 MeV. This system has been 
extensively studied over the years and a number of serious problems has re- 
mained unsolved: The explanation of the anomalous large angle scattering 
data; the out-of-phase problem between theoretical predictions and experi- 
mental data; the reproduction of the oscillatory structure near the Coulomb 
barrier; the consistent description of angular distributions together with the 
excitation functions data are just some of these problems. We propose the 
use of a modified potential method to explain these problems over this wide 
energy range. This new method consistently improves the agreement with 
the experimental data and achieves a major improvement on all the previous 
Optical model calculations for this system. 
Keywords: 16 0+ 28 Si Reaction, Optical model, elastic and inelastic scattering, anoma- 
lous large angle scattering (ALAS), excitation function. 
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I. INTRODUCTION 



The elastic and inelastic scattering between heavy ions have been amongst the main 
sources of information about complex nuclei over the last 40 years. [1-3]. There have been 
numerous experimental investigations for the systems with a combined mass number of 
A P + A T < 60, and these investigations have displayed a common unexpected feature near 
#cm=180° for the elastic and inelastic scattering cross-sections. 

The physical origin of the observed structure is not yet fully understood [2,4-8] and 
presents a challenge to different approaches that have been proposed to explain it. These 
approaches range from the occurrence of possibly overlapping shape resonances [9] and the 
scattering from surface-transparent Optical potentials [10] to more exotic effects like explicit 
parity-dependence of the ion- ion potential [11,12]. The first approach that has been more 
popular with researchers so far attempts to describe the data by invoking the properties of 
the average Optical potential [10,13-19]. On the other hand, in the second approach, the 
structure in the excitation function is associated with isolated or nearly isolated partial- wave 
resonances superimposed on the scattering properties of a standard Optical potential [1,2]. 
At present, none of these approaches provides a consistent explanation for all the existing 
data for this system. 

Consequently, the following problems continue to exist for this reaction [5,8,20,21]: (1) 
The explanation of anomalous large angle scattering data; (2) the reproduction of the oscil- 
latory structure near the Coulomb barrier; (3) the out-of-phase problem between theoretical 
predictions and experimental data; (4) the consistent description of angular distributions 
together with the excitation functions data. 

Within the framework of the coupled-channels method, we have analyzed such light 
heavy-ion reactions that pose the above-mentioned problems [5,8,20,21]. The new approach 
that we proposed within the coupled-channels method has successfully explained the ex- 
perimental data over wide energy ranges for different systems such as 12 C+ 12 C, 16 0+ 28 Si, 
16 0+ 24 Mg and 12 C+ 24 Mg. One feature we observed in these analyzes was that these reac- 
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tions were extremely sensitive to the shape of the nuclear potential in the surface region. In 
this paper, by taking this feature into account, we consider an extensive investigation of the 
elastic scattering of this system at numerous energies by using a modified potential. Similar 
to the previously conducted coupled-channels analyzes for such systems [5,8,20,21], this new 
technique modifies the shape of the potential at the surface. Thus, we aim to address the 
above-mentioned problems within the framework of the Optical model and to obtain results 
that are as good as the results of the coupled-channels method. Accordingly, we analyze the 
experimental data from 29.0 to 142.5 MeV in the laboratory system over the whole angular 
range up to 180°. The 180° elastic scattering excitation function has also been studied over 
this energy range. Extensively modified version of the code CHUCK [22] has been used for 
the all calculations. 

In the next section, we introduce our Optical model and potential parameters to explain 
the observed experimental data. Then, we show the results of these analyzes in section III 
from E iafe =29.0 MeV to 142.5 MeV. Our conclusion is given in Section IV. 



The standard Optical model with folding model potentials or with similar phenomenolog- 
ical potentials such as the square of the Woods-Saxon has failed to describe certain aspects 
of the experimental data. Therefore, similar to Mackintosh, Kobos and Satchler's work 
[17-19], in the present calculations, our total real potential consists of the nuclear potential, 
VN UC i ear , with two small additional potentials, [U(r) = Ui(r) + U 2 (r)]: 



The nuclear potential is assumed to have the square of a Woods-Saxon shape and the 
parameters are fixed to reproduce the folding model potential of Mackintosh, Kobos and 
Satchler [17-19]: 



II. THE MODEL 



V t otal(r) = V Nudear (r) +U(r) + V C oulomb(r) + V Ce ntrifugal(r) 



(1) 



Real potential 




-v 



(2) 



(1 + exp{r — R)/a) 2 
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where ^=761.5 MeV and R=r (A a 1 / 3 +A A 1 / 3 ) with r =0.75 fm and a=1.425 fm. The pa- 
rameters of the nuclear potential are fixed as a function of energy and kept constant in the 
present calculations although small changes were observed to improve the quality of the fits. 

The Coulomb potential [23] due to a charge Z a e interacting with a charge Z A e, distributed 
uniformly over a sphere of radius R c , is also added. 

V Co uiomb{r) = _L Za ^ e r>R c (3) 

47re r 

1 Z a Z A e\ r\ 

"(3 - -^), r < R c (4) 



47re 2R C v R 2 C 

where R c =7.78 fm is the Coulomb radius, and Z a and Z A denote the charges of the projectile 
a and the target nuclei A respectively. 

The sum of the nuclear, Coulomb and the centrifugal potentials is shown in figure 1 for 
various values of the orbital angular momentum. The superposition of the attractive and 
repulsive potentials results in the formation of a potential pocket. The width and depth of 
the pocket depend on the orbital angular momentum quantum number for a given nuclear 
potential. This pocket is very important for the interference of the barrier and internal 
waves, which creates the oscillatory structure observed in the cross-section. 

The significance of the two small additional potentials should be emphasized here, since 
without them, it would be impossible to fit the experimental data in the elastic scattering 
calculations. These two small additional potentials are the derivatives of the Woods-Saxon 
shape and the parameters are shown in table I where 

JJ ( r ) = 4[/ 1 a 1 d/(r, R u oi) /dr + AU 2 a 2 df(r, R 2 , a 2 ) /dr (5) 

f(r, Ri, cii) = — -4 — — - (6) 

(1 + exp{{r - Rij/ai)) 

As shown in the insert of figure 2, they create two minima in the nuclear potential 
between ~5.5 fm and ~9.0 fm. The effect of these potentials can be understood in terms of 
the interference between the internal and barrier waves that correspond to a decomposition 
[24,25] of the scattering amplitude into two components, the inner and external waves. The 



inner wave comes from the reflection at the inner face of the total real potential pocket and 
the external wave comes from the reflection at the outer barrier (see figure 1 for the pocket 
in the total real potential). The presence of the two small potentials affects the phases and 
magnitudes of these internal and external components. We observe from the parameters in 
table I that the two small additional potentials are not strong enough to produce pockets 
in the total real potential although they have a very significant effect on the scattering. 
This is demonstrated in figure 4 at 41.17 MeV. In that figure, (d) shows the cross-sections 
obtained when both potentials are omitted. As a result, the calculated oscillations are out- 
of-phase with the measured ones even at the intermediate angles. Without these potentials, 
we were unable to refit the data by merely varying the parameters of the real and imaginary 
potentials. Thus, the oscillatory structure could not be reproduced correctly. 

The imaginary part of the potential was taken as the sum of a Woods-Saxon volume and 
the surface potential [19]: 

W(r) = -W v f(r, R v , a v ) + 4W s a s df(r, R s , a s )/dr (7) 

fMa)= (l + exp((r-R)/a)) (§) 

with WV=59.9 MeV, a v =0.127 fm and W s =50.0 MeV, a 5 =0.250 fm. These parameters 
were also fixed in the calculations and only their radii were increased linearly with increasing 
energy according to the following formulae. 

R v = 0.06084£ CM - 0.544 (9) 

R s = 0.2406£ CM - 2.191 (10) 

The imaginary potentials are shown in figure 3 for E Lab — 41.17 MeV. 

Moreover, the relative significance of the volume and surface components of the imaginary 
potential was evaluated for all the energies. For higher energies, omitting the volume term 
predominantly affected the amplitude of the cross-section at large angles. However, this 
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effect was small and negligible at lower energies. Omitting the surface term increased the 
cross-sections at large angles which were as much as two orders of magnitude. It was observed 
that this term had a significant effect at all the considered energies. 

III. THE RESULTS 

Using the above-described Optical model with fixed real and linearly increasing imaginary 
potential parameters, we have analyzed the experimental data from 29.0 MeV to 142.5 MeV 
in the laboratory system. The comparisons between experimental data and the Optical 
model fits are shown in figures 5, 6 and 7. The x 2 values for some of the energies studied 
are also shown in table III. As it can be seen from these figures and the x 2 table, we 
have obtained excellent agreement with the experimental data over the whole energy range 
considered. 

We have also analyzed the averaged value of the excitation function over the angular 
range 180° ±5° with these potentials. The radii of the imaginary potential increased linearly 
with energy up to 54.0 MeV according to equations (9) and (10). Beyond this energy, we do 
not have the available experimental data for the large angles. Therefore, as shown in table II, 
the radii were kept constant at that energy and were not changed for the higher energies. The 
result is compared with the measured values in figure 8. The overall agreement is quite good 
up to about 50 MeV where experimental data are available at large angles. In particular, the 
maxima and minima in the excitation function are obtained correctly. However, it should 
be emphasized that the calculated excitation function is much more sensitive to the details 
of the Optical potential than are the fits to the angular distributions. For example, omitting 
two small potentials completely destroys the agreement with the excitation function. 

We were able to obtain almost the same agreement with the experimental data above 
50 MeV regardless of these two small potentials. This shows that for higher energies, the 
barrier /internal wave decomposition is not very important and thus we may remove these 
two small potentials. By using the same nuclear potential and the same type of imaginary 
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potentials, but by readjusting the depth and diffuseness of its surface component, we were 
able to fit the experimental data from E Lab =50 MeV to 142.5 MeV. These results are shown 
in figure 7 and the relevant parameters are shown in table II. 

IV. SUMMARY 

We have shown a consistent description of the elastic scattering of the 16 0+ 28 Si system 
from 29.0 MeV to 142.5 MeV in the laboratory system by using the Optical model calcu- 
lations. In the introduction, we presented the problems that this reaction manifests. We 
attempted to find a consistent solution to these problems. However, within the standard 
Optical model calculations, we failed, as others did, to describe certain aspects of the data. 

Our previous elastic and inelastic analyzes with light-heavy ion reactions using the 
coupled-channels formalism had already shown that these reactions are extremely sensi- 
tive to the shape of the potential in the surface region. By taking this feature into account, 
we have used two small additional potentials that modify the shape of the potential in the 
surface region. As a result, we have obtained excellent agreement with the experimental 
data over a wide energy range similar to that of the coupled-channels method. The compar- 
ison of the results indicates that a global solution to the problems relating to the scattering 
observables of this reaction over a wide energy range has been provided by this method. 
Finally, it should be noted that although these additional two small potentials we used are 
very small and do not create a pocket in the total nuclear potential, they are very effective 
for the interference of the barrier and internal waves, which creates the oscillatory structure 
observed in the cross-section. Further work in order to derive this term from a microscopic 
viewpoint is still under-progress. 
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TABLES 



Ui 


Ri 


ai 


u 2 


i?2 a 2 


(MeV) 


(fm) 


(fm) 


(MeV) 


(fm) (fm) 


10.14 


6.057 


0.177 


3.17 


7.358 0.371 




TABLE I. The parameters 


of the two small additional potentials. 


W v 


R v 


ay 


w D 


Rd a D 


(MeV) 


(fm) 


(fm) 


(MeV) 


(fm) (fm) 


59.90 


1.3925 


0.127 


31.25 


5.4697 0.550 



TABLE II. The parameters of the potentials required to fit the higher energy data. These 
are the values of the imaginary potential at £x a fc=50 MeV and they were fixed to reproduce the 
structure for forward angles at high energies for 16 0+ 28 Si. 
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Laboratory Energy 



Optical Model 



29.34 


1.2 


29.92 


1.4 


30.70 


1.9 


31.63 


1.5 


32.75 


0.3 


33.17 


0.8 


33.89 


0.9 


35.04 


2.0 


35.69 


8.0 


38.20 


7.9 


41.17 


16.5 



TABLE III. x 2 values for some of the energies studied. 
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FIGURES 
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FIG. 1. Interaction potential between O and Si is plotted against the separation R for 
various values of the orbital angular momentum quantum number, I. The parameters are given in 
the text. 
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FIG. 2. The shapes of two small additional potentials U\ and XJi are displayed by dotted and 
solid lines respectively. The inserted figure shows their effects on the nuclear potential with a 
long-dashed line. 
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FIG. 3. The volume and the surface components of the imaginary potential at £^{,=41.17 
MeV for the 16 0+ 28 Si system. 
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FIG. 4. The effect of the two small additional potentials for the £x a b=41.17 MeV. (a) reflects 
that the best fit is obtained with the inclusion of both potentials, (b) displays the situation without 
the inclusion of U2 and (c), without the inclusion of U\. Finally, (d) displays the outcome when 
both potentials are ignored. 
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FIG. 5. The elastic scattering results obtained by using the single-channel Optical model cal- 
culations for the 16 0+ 28 Si system. 
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FIG. 6. Elastic scattering results obtained by using the single-channel Optical model calcula- 
tions for the 16 0+ 28 Si system (continued from figure 5). 
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FIG. 7. Results of the elastic scattering calculation for forward angles at higher energies by 
using the single-channel Optical model calculations. 



18 



10° 



10" 1 

c? 

-o 
D 

10" 2 
10" 3 

20.0 30.0 40.0 50.0 60.0 70.0 80.0 90.0 
Laboratory Energy (MeV) 

FIG. 8. Comparison of the calculated 180° elastic scattering excitation function with the mea- 
sured experimental data. The theoretical calculation is averaged over the last 5°. 
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